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Abstract. A new class of geometric query problems are studied in this paper. We are required to 
preprocess a set of geometric objects P in the plane, so that for any arbitrary query point q, the largest 
circle that contains q but does not contain any member of P, can be reported efficiently. The geometric 
sets that we consider are point sets and boundaries of simple polygons. 

Keywords: Largest empty disk, query answering, medial axis, computational geometry 
1 Introduction 

Largest empty space recognition is a classical problem in computational geometry, and has 
applications in several disciplines like database management, operations research, wireless 
sensor network, VLSI, to name a few. Here the problem is to identify an empty space of a 
desired shape and of maximum size in a region containing a set of obstacles. Given a set 
P of n points in M^, an empty circle, is a circle that does not contain any member of P. 
An empty circle is said to be a maximal empty circle (MEC) if it is not fully contained in 
any other empty circle. Among the MECs, the one having the maximum radius is the largest 
empty circle. The largest empty circle among a point set P can easily be located by using 
the Voronoi diagram of P in 0{nlogn) time |i32j . 

Although a lot of study has been made on the empty space recognition problem, surprisingly, 
the query version of the problem has not received much attention. The problem of finding the 
largest empty circle centered on a given query line segment has been considered in [3]. The 
preprocessing time, space and query time complexities of the algorithm in [3] are 0{n^ logn), 
O(n^) and O(logn), respectively. In practical applications, one may need to locate the largest 
empty circle in a desired location. For example, in the VLSI physical design, one may need 
to place a large circuit component in the vicinity of some already placed components. Such 
problems arise in mining large data sets as well, where one of the objectives is to search for 
empty spaces in data sets j2^. In [12], Edmonds et al. formalized the problem of finding 
large empty spaces in geometric data sets. In particular, they studied the problem of finding 
large empty rectangles in data sets. 
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An important problem in this context is the circular separability problem. Two planar sets Pi 
and P2 are circularly separable if there is a circle that encloses Pi but excludes P2. O'Rourke 
et al. [27] showed that the decision version of the circularly separability of two sets can 
be solved in 0{n) time using linear programming. Furthermore, they show that a smallest 
separating circle can be found in 0{n) time while the computation of the largest separating 
circle needs O(nlogn) time. Detailed study on circular separability problem can be found 
in [7f5|13ll27j . Boissonnat et al. [5] proposed a linear-time algorithm for solving the decision 
version of circular separability problem where the sets Pi and P2 are simple polygons, and 
the algorithm outputs the smallest separating circle. They also consider the query version 
of this problem where the objective is to preprocess a convex polygon P such that given a 
query point q and a query line £, report the largest circle inside P that contains q and does 
not intersect ^. The preprocessing time and space complexities of their proposed algorithm 
are both 0(r;,logn), and the query can be answered in 0(log?T,) time. They also showed that 
a convex polygon P can be preprocessed in 0{n) time and space such that for a query set S 
of k points, the largest circle inside P that encloses S can be computed in 0(fclog?T.) time. 

In addition to empty circles, empty rectangles have also been studied. We introduced the 
query version of the maximal empty rectangle in [1] . The problem entails preprocessing a set 
of n points such that, given a query point g, the largest empty rectangle containing q can be 
reported efficiently. We gave a solution with query time OiXogn) with preprocessing time and 
space being 0(n^ log n) and O(n^), respectively. Recently, Kaplan et al. [19jJ improved the 
preprocessing time and space complexities to 0{na{n)\og^ n) and 0(na(ra) log^ n), respec- 
tively, while the query time has increased to O(log^n). Here a{n) is the inverse Ackermann 
function. 

1.1 Our Results 

In this paper, we study the query versions of the maximum empty circle problem (QMEC). 
The following variations are considered. 

— Given a simple polygon P, preprocess it such that given a query point g, the largest circle 
inside P that contains the query point q can be identified efficiently. 

— Given a set of points P, preprocess it such that given a query point g, the largest circle 
that does not contain any member of P, but contains the query point q can be identified 
efficiently. 

Our results are summarized in Table [H 

We believe that our work will motivate the study of new types of geometric query problems 
and may lead to a very active research area. The main theme of our work is to achieve 
subquadratic preprocessing time and space, while ensuring polylogarithmic query times. 
The results in this paper, improve upon the results in our previous work [T|. Very recently, 
Kaplan and Sharir [20j provided a solution to the QMEC problem for point sets that only 
requires O(nlog^n) time and 0{n\ogn) space for preprocessing. Their query times, however, 
are O(log^n). 
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Table 1. Complexity results for different variations of largest empty space query problem 



Geometric set 


Preprocessing 
time 


Space 


Query time 


Simple Polygon 


0(n log^ n) 


0(n log ri) 


O(logn) 


Point Set 


0(^3/2 log^ n) 


0(^3/^ logn) 


(log n log \ogn) 


Point Set 


0{n^/^ \ogn) 




O(logn) 



1.2 Organization of the paper 

In Section [2| as a preliminary requisite, we describe a way to answer QMEC query for the case 
of convex polygons. The same bounds have been achieved by Boissonnat et al. [8J, but our 
solution is slightly different and serves as the basis for solving the QMEC problem on simple 
polygons. In Section [3} we present the QMEC problem for simple polygons with n vertices. 
The preprocessing time and space complexities are 0{n\og^ n) and 0{n\ogn) respectively, 
and the query answering time is 0(log?T,). In Section |4| we consider the same problem on a 
set P of n points in M^. We present two algorithms (cf. Table [T]). Our first algorithm uses 
the concept of planar separators [22] on the underlying planar graph corresponding to the 
Voronoi diagram of P. It solves the QMEC problem on P with ©(n^/^log n) preprocessing 
time and 0{n^^^\ogn) space. Here, the queries can be answered in O (log n log log n) time. 



Our second algorithm (cf. Section 4.4) uses the r-partitioning [T3] of planar graphs. With a 
suitable choice of r, the query time is only O(logn), an improvement over our first algorithm. 
However, the preprocessing time and space increase to 0{n^^^ logn) and 0{n^^^), respectively. 



2 Preliminaries: QMEC problem for a convex polygon 

Let P be a convex polygon and {pi,P2, . . . ,Pn} be its vertices in counter-clockwise order. Our 
objective in this section is to preprocess P such that given an arbitrary query point q E P, 
the largest circle Cg containing q inside the polygon P can be reported efficiently. Note that, 
the locus of the centers of all the maximally empty circles (MECs) inside P is defined to be 
the medial axis M of P. Let c be the center of the largest MEC inside P (see Figure [l]^a))|^ 
The medial axis of a convex polygon consists of straight line segments and can be viewed 
as a tree rooted at c [11]. To avoid confusion with the vertices of the polygon, we call the 
vertices of M as nodes. Note that, the leaf- nodes of M are the vertices of P. Let us denote 
an MEC of P centered at a point x G M as MEC^, and let be the area of MEC^;. 

In [8] , a planar map of circular arcs is constructed by drawing the MEC at each node of M in 
0{n) time and space. The problem of finding Cg reduces to the point location problem in the 
associated planar map. These point location queries can be answered in O(logn) time. We 
propose an alternative solution (with the same complexity results as in [S]) because our new 
technique plays a basic role in solving the problem when P is a simple polygon (cf. Section 

^ There can be infinitely many MECs of largest radius, in which case we pick c to be the center of one such MEC. 
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|3.2[ ). We use the fact that the medial axis M is a tree, and then use the level-ancestor queries 
[6] on M. 




Fig. 1. (a) Illustration of Lemma [Tl and (b) partition of P. 



Lemma 1. fBj As the point x moves from the center c of the largest MEC along the medial 
axis towards any vertex pi G P (leaf node of M), decreases monotonically. 

Lemma 2. The polygon P can be preprocessed in 0{n) time such that given any arbitrary 
query point q, a point x on M such that MEC-r contains q can be reported in 0{logn) time. 

Proof. The medial axis M subdivides P into n convex faces such that each face Pi consists 
of an edge PiPi+i from P and a convex chain of edges from M connecting pi to pi+i (see 
Figure [l](b)). In the preprocessing phase, we perform the following steps. 

1. Compute the medial axis M of P, which is a tree rooted at c. 

2. Compute the subdivision in 0{n) time. For this we will need M, which can be computed 
in linear time [TT] . 

3. Store the chain of edges associated with each face in an array so that it is amenable to 
binary searching. 

4. Finally, the subdivision can be preprocessed in 0{n) time so that the face containing a 
query point q can be located in O(logn) time [2T]. 

In the query phase, we perform the following steps. 

1. We find the face Pi that contains q in O(logn) time. 

2. Recall that exactly one edge PiPi+i of P will be an edge in Pj. Consider the line i through 
q that is perpendicular to the edge PiPi+i- It will intersect an edge in M that is also an 
edge bounding P^; we report that intersection point as x. Note that x can be computed 
in O(logn) time via binary searching over the chain of medial axis edges bounding Pj. 
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We need to prove that MECj; indeed encloses q. Firstly, note that (. will intersect the edge 
PiPi+i internally at a point t because (i) Pi is convex and (ii) the two internal angles in Pi 
at Pi and pi+i are both acute. Secondly, note that any MEC that goes through t must be 
tangential to PiPi+i, thereby making it unique and centred on i; more precisely, the MEC that 
goes through t must be centred at x. Finally, from the construction, it is clear that q lies on 
the diameter of MEC^, thus proving that that MEC^,. encloses q. □ 

Now we will describe how to solve the QMEC problem for a convex polygon. Given a query 
point q we find (using Lemma [2]) the point x on M such that MEC^,. encloses q. 

Observe (informally for now) that, for any fixed point q inside P, the MECs that encloses q are 
centered on a connected subtree iVf of the medial axis M. This observation is formally proved 
in Lemma |3] in the more general setting of simple polygons. Coupling this observation with 
Lemma [l| we can conclude that Cq is the point on M'^ closest to the root c of M. Therefore, 
we can locate Cq by performing a binary search on the path x ~ c. We find two consecutive 
nodes v and its parent v' on the path x ~ c such that MEC^, encloses q, but MEC„/ does not. 
Since the path lies on a tree representing the medial axis M, we can use level-ancestor queries 
[6] for this purpose. After computing v and v' , the exact location of Cq can be computed in 
0(1) time. Thus, we have the following theorem: 

Theorem 1. A convex polygon P with n vertices can be preprocessed in 0{n) time and space 
such that, given any arbitrary query point q E P, the largest circle containing q inside P can 
be reported in O(logn) time. 



3 QMEC problem for simple polygons 

Let P be a simple polygon on n vertices. Recall that the medial axis M of P is defined to be 
the locus of the centers of all circles inside P that touch the boundary of P in two or more 
points (see, e.g., [H]). While the medial axis of a convex polygon consists only of straight 
line segments, the medial axis of a simple polygon may additionally contain parabolic arcs 
[28]. 

Our approach for solving the QMEC problem uses the fact that M is a geometric tree. Its leaf 
nodes correspond to the vertices of P, and the internal nodes correspond to the points on M 
such that the MECs centered at each of those points touch three or more distinct points on 
the boundary of P. We denote the set of internal nodes of M as N. An edge in M is a path 
between two nodes that does not contain any other node in its interior. Note that a single 
edge consists of one or more line segments or parabolic arcs. 

For any point x G M, we denote the maximal empty circle centered at a; in P by MEC^-. A 
point X G M, that is not a leaf, is said to be a valley (resp., peak) if for a positive 5 — )■ 0, 
the MECs centered at points in M within distance 5 from x are at least as large as (resp., no 
larger than) MEC^^ and at least one such MEC is strictly larger (resp., smaller) than MEC^. Note 
that a pair of parallel edges in P may induce a pair of peaks or a pair of valleys. In such 
cases, we only pick one representative peak or valley and discard the other. We use ^ and 
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O to denote the set of valleys and peaks, respectively. For any x G it is easy to observe 
that ViECx touches P in exactly two points diametrically opposed to each other. Otherwise, 
we can move along a direction to get smaller MECs. As a consequence, a valley can only be 
in the interior of an edge. Therefore, <P fl N = 0. On the other hand, 6* C N. 

We define a mountain to be a maximal subtree of M that does not contain any valley point 
(except as its leaves). We partition M by cutting the tree at all the valley points, and generate 
a set of mountains M = {Mi, M2, . . . , M|m|} (See Figure |2|^a)). 

Observation 1 

(i) Each valley point is the common leaf of exactly two mountains, 
(a) Each mountain has exactly one peak. 

(Hi) If a point x moves from a valley point of a mountain towards its peak, the size of KEC^ 
increases monotonically. 

Proof. Suppose f is a valley point. We have noted earlier that v can only be in the interior 
of an edge. Therefore, v must be a common leaf between exactly two mountains. 

For part (ii), assume that there are two peaks pi and p2 in a mountain. Consider the path 
in M from pi to p2 (denoted pi ~ ^2)- Consider the point x* = argmin^-gp^^pj MECa;. One can 
observe that x* will be a valley, implying that pi and p2 cannot be in the same mountain. 

Consider a point x that moves from a valley of a mountain towards its peak. If the MECs 
don't increase monotonically, it is easy to see that a valley point will be encountered. This 
implies that x has moved into another mountain. □ 



Fig. 2. (a) Partitioning the medial axis M into mountains, and (b) The subtree M'' for a query point q 

At each valley point x of M, consider the chord in P connecting the two points at which 
MECa; touches P. These chords induced by each x E 'P will partition P into a set of sub- 
polygons {Pi, P2, • • • , -P|M|} of cardinality equaling the total number of mountains because 
this partitioning ensures that the portion of M contained in each of these sub-polygons is a 
mountain containing a single peak. 

Given a (query) point g G P, let M'^ C M denote the locus of the centers of all possible 
maximal empty circles in P that enclose q (see Figure |2](b)). The following structural lemma 
plays a crucial role in designing our algorithm. 
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Lemma 3. M'^ is a (connected) subtree of M. 

Proof. For the sake of contradiction, let us assume that M'^ is disconnected. Then, there are 
at least two distinct points ti and t2 on the medial axis such that (a) MECfj and MECtj contain 
q, but (b) MECt does not contain q for some t on the path along the medial axis connecting 
ti and ^2- 

Without loss of generality, we assume that such a t is not a node in M. Therefore, MECf 
touches the simple polygon P at exactly two points, a and b. The chord [a, b] partitions P 
into two polygons Pie ft and Pnght to the left and right of [a, b] respectively (see Figure |3](a)). 
Note that, t also partitions the medial axis into two subtrees, M/g/j and Mj-ight, such that 
ti G Mieft and ^2 ^ Mright- For the rest of the proof, we use MECj and Pieft to denote the 
region enclosed by them. We now claim that 

mCt, C Pieft U MEQ = P \ {Pright \ MEQ). (1) 




(a) (6) 

Fig. 3. Proof of Lemma [s] (a) the general case, and (b) a special case where the points a and b are concave vertices 
of P 



To show that Equation [T] holds, we consider the following cases. 

Case: MECt^ touches P at both a and b: This case is illustrated in Figure [sj^b). Here, a 
and b must be concave vertices that induce a straight line segment edge in the medial 
axis. Both t and ti are on that edge; in particular, ti will be to the left of t. It is now 
easy to infer that Equation [T] holds. 

Case: MECt^ touches at most one of {a, b}: Let MECtj touch the other point d ^ {a, b} on 
the boundary of P. Clearly, d G Pie/t \MEC(. If we assume that MECj^ also passes through a 
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point d! G Pright\^^'^ti then it is impossible to construct MECj^ without properly enclosing 
some point outside P. Therefore, Equation [T] holds. 

By symmetry, we can also say that 

Pright UMECt = P\(Pie/AMECt). (2) 
Taking the intersection of Equations [l] and [2| we get 

MECt, n MECt, C {Pieft U MECt) H {Pright U MECt) 

= (Pje/t n Pright) U MECt = MECt (since (P^^/t n Pr^ght) C MECf 

This contradicts our assumption that q falls in MECt, and MECt,, but not in MECt- □ 

Corollary 1. Let v E M be such that MEC^ does not contain q. Then is contained entirely 
in one of the subtrees of M obtained by deleting v from M . 

Corollary 2. Consider two points u, f G M, such that MECt, overlaps with MEC^- Let a be a 
point in P such that a G MEC^, fl MEC„. Then a will lie in all MECs centered along the path 
from V to u in M. 

Proof. Since MECt, overlaps with MEC„ and a G MEC„ fl MEC^, both v and u are points on M". 
The result now follows immediately from Lemma |3} □ 

Before we delve into solving QMEC, in the next three subsections, we define three data struc- 
tures that we use as building blocks. 

3.1 PLiCA: Point location in circular arrangement 

The problem is to preprocess a set C = {Ci,C2, ■ ■ ■ ,Cn} of circles of arbitrary radii, so 
that for any query point q in the plane, we need to quickly report if there exists a circle 
Cj G C such that g G Cj. This can be achieved by using the concept of Voronoi diagrams in 
Laguerre geometry of circles in C [16]. Each cell of the Voronoi diagram is a convex polygon 
and is associated with a circle in C. The membership query is answered by performing a point 
location in the associated planar subdivision. The preprocessing time and space complexities 
are 0(?T,logn) and 0{n) respectively, and the queries can be answered in 0(logr;,) time. 

3.2 QiM: Query-in-Mountain 

Given a mountain Mj G M, we must preprocess it such that given a query point q inside P 
such that Mj n Af? 7^ (and a point x G n M'^), our task is to report the largest MEC 
centered at a point on M, that contains q. Note that if the center moves from x towards the 
peak of Mi, the size of the MEC increases monotonically. Thus, we can apply the algorithm 
proposed in Section |2] for the convex polygon case to identify the largest MEC containing q, 
and centered on Mj fl M'^. The preprocessing time and space complexities for the mountain 
Mi are both 0(|Pj|), and the query time is 0(log|Pj|), where \Pi\ denotes the number of 
edges in the sub-polygon Pi that induces Mj. Since the set of mountains and sub-polygons 
are partitions of the medial axis M and the polygon P, respectively, all the mountains can 
be preprocessed for the QiM queries in 0{n) time. 
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3.3 QiC: Query-in-Circle (Problem Definition and Bounds) 



The QiC is a simplification of the QMEC problem in which, in addition to P and its medial axis 
M, a node w of M is specified as part of the input for preprocessing. We are promised that 
the query point q will lie inside MEC^. As in QMEC, we are to report the largest MEC Cg that 



contains q. We defer the details of our solution for the QiC problem to Section 3.6, where we 
prove the following theorem. 

Theorem 2. There exists a solution for QiC that takes 0{nlogn) time and 0{n) space for 
preprocessing. Queries can he answered in 0{\ogn) time. 

To solve QMEC, we employ a divide-and-conquer strategy that divides the medial axis into 



smaller pieces. On these smaller pieces, we employ QiC. We remark in Section [3^ how the 
solution to the QiC problem on the entire medial axis can be adapted to restricted portions 
of the medial axis. For now, we note that the preprocessing time and space of QiC scale with 
the size of the portion of the medial axis that is preprocessed. On a portion M* of the medial 
axis, the preprocessing time and space are 0{n* \ogn*) time and 0(n*), respectively, where 
n* is the number of edges of P that induc^the edges in M* (cf. Corollary [4]) . 

Algorithm 1 Preprocessing for QMEC on a Simple Polygon P 
Require: A simple polygon P. 

Compute the medial axis M of P. 

Construct a PLiCA data structure on MECs centered on nodes of M. 
Construct a secondary PliCA data structure on the MECs centered on valley points of M. 



4: Construct a list (Mi, M2, . . .) of mountains and preprocess each mountain for QiM. 
5: Partition P into sub-polygons Pi, P2, ■ ■ ■ such that each Pj is associated with its cor- 
responding Mi. (Recall that this can be performed by cutting along diameters of MECs 
centered on valley points.) 
6: Preprocess P and its sub-polygons (in 0{n) time and space) such that, given a query 
point g, the sub-polygon that contains q can be reported efficiently (in O(logn) time). 
Call this data structure D. 
T <«— Decompose(M). 
for each node t & T do 

Let M* C M be the subtree associated with t. 
10: Let f* be the centroid of M* (cf. Lemma |4]). 

11: Preprocess MEC^t for QiC with the additional promise that the largest empty circle Cq 
that contains query point q is centered on M*. Associate this QiC data structure with 
t. 

12: end for 



We say that an edge ep in P induces an edge cm in M if for some point x in the interior of cm, KECx touches ep. 
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Algorithm 2 Decompose(M) 



Require: A tree M with n nodes. 

Ensure: A divide-and-conquer tree T that decomposes M. 
1: if n = 1 then 

2: return a tree with the single node. 
3: end if 

4: Find the centroid (cf. Lemma [4| of M that will decompose M into subtrees M^, M^, . . . 

5: Create a node node with a list child of child pointers. 

6: Associate M with node. 

7: for each subtree do 

8: node . child [i] ^ Decompose (M*). 

9: end for 



Algorithm 3 Query phase of QMEC on a simple polygon P with query point q 
Require: Query point q and all the data structures created in the preprocessing phase. 
{We can use the PLiCA data structure for the following condition.} 

1: if q falls inside some MEC centered on a node v oi M then 

2: {We are in the affirmative case.} 

3: Find the node t in T whose centroid is v. 

{In the next step, v* is the centroid of the subtree of M associated with t* .} 

4: Find (in T) the farthest ancestor t* of t such that MEC^* contains q. 

5: Cq ^ circle returned by querying the QiC data structure associated with v* . 

6: Return Cq. 

7: else 

8: {We are in the negative case.} 

9: if q falls in some MEC centered on a valley point p then 
10: {p is a valley point connecting exactly two mountains A and B.} 
11: Cq is the larger of the MECs returned by querying the QiM data structures associated 

with A and B. 
12: Return Cq. 
13: else 

14: Use data structure D to find the sub-polygon P'^ that contains q. 

15: Cq i- circle returned by querying the QiM data structure associated with P'^. 

16: Return Cq 

17: end if 

18: end if 
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3.4 Preprocessing for the QMEC problem 



Algorithm [T] outlines the steps in the preprocessing phase. The first [6] steps are straight- 
forward. Before we explain the subsequent steps, we state (for the sake of completeness) a 
well-known lemma. 

Lemma 4. Every tree M with n nodes has at least one node v whose removal splits the 
tree into subtrees with at most [|] nodes. The node v is called the centroid of M . 

In line number [T] we call Algorithm |2] (recursively) to build a centroid decomposition tree T. 
We partition using the centroid (cf. Lemma |4]) in order to ensure that T is balanced. 

The centroid decomposition is constructed in anticipation of the query phase. Suppose q is 
a query point. If q lies in MEC„, where v is the centroid associated with the root of T. Then, 



we can use the QiC attached to the root (in line number 11). If, on the contrary, q ^ MEC^, 
then from Corollary [T] we know that only one of the subtrees rooted at v will contain M'^, 
thereby allowing us to recurse into that subtree until we find the centroid whose MEC encloses 
q. To facilitate this recursion, we must provide a way to find the correct subtree to recurse 
into. For this, we consider the geometry of the polygon P. Let MEC„ touch the polygon P 
at fc (> 3) points. Consider the partitioning of P into k sub-polygons, apart from the one 
containing v, by inserting chords as shown in Figure |4} These k sub-polygons correspond to 
the k subtrees of M obtained by removing f . It is easy to see now that a point location data 
structure will suffice. In the query phase, we can simply find the sub-polygon that contains 
q and recurse into the corresponding subtree. 




Fig. 4. The divide and conquer search structure 



In lines [s] to [12} for each node t of T we associate an appropriate subtree M* of M along with 
the centroid f * of M*. Additionally, we will construct a QiC data structure associated with t 
with the additional promise that the largest empty circle Cg that contains query point q is 
centered on M*. 
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Lemma 5. The time and space required for preprocessing P are 0{n\o^ n) and 0{n\ogn), 
respectively. 

Proof. The medial axis M of a simple polygon can be computed in 0{n) time [11]. Once 
we have M, the partition of M into mountains, and the associated partitioning of P can be 
done in 0{n) time. The data structure for the planar point location can easily be obtained 
in 0(n log n) time and 0{n) space. The PLiCA data structure for all the MECs centered at the 
nodes of M requires 0(?T,log?7,) time and 0{n) space. 

Consider a level ^ in the tree T. Each node in level i implements the QiC data structure 
on a portion of the medial axis that is disjoint from the portion addressed by other QiC 
implementations in the same level. Therefore, the preprocessing times and spaces of all QiCs 
at any level £ is 0{n\ogn) and 0{n) respectively. Since there are O(logn) levels in T, the 
total preprocessing time and space for all QiCs is 0{n log^ n) and 0{n log ra), respectively. □ 



3.5 QMEC query 



As discussed in Algorithm [3} in the query phase with a point g, we first test whether q lies 
in the MEC centered at any node v of the medial axis. This can be performed using the PLiCA 



data structure (see Subsection 3.1 ) built over the set of MECs centered at all the nodes in M. 
We now need to consider two cases: 

Affirmative case: There exists a node f in M such that MEC^, contains q. From t ^ T 
corresponding to v we move upward in the centroid tree T following the parent pointers to 
identify a node t* E T at the highest level such that the MEC centered at v*, the medial axis 
node of the subtree associated with t*, contains q. Our choice of v* coupled with Corollary 
[l] imply the following lemma. 

Lemma 6. Let M* be the subtree of M associated with t* and let v* be the centroid of 
M*. Then, W C M* and MEC„* contains q. 

Lemma [6] ensures all the prerequisites for QiC, so we can query the QiC data structure 
associated with node t* and correctly obtain Cq. 

Negative case: There exists no node f in M such that MEC„ contains q. In this case, M'^ 
cannot span more than two mountains as otherwise M'^ must include a node in M. If q falls 
in an MEC centered at a valley point p, we query the QiM data structure associated with the 
two polygons connected by p. Otherwise, the center of Cq lies in only one mountain, the 
mountain that contains q. We identify the sub-polygon Pj in the planar subdivision that 
contains q using the data structure D (see line number |6] of Algorithm [T]) . Finally, we can 
compute Cq by performing the QiM query in Mj, the mountain associated with Pj. 

Theorem 3. Given a simple polygon P, we can preprocess it in 0{n log^ n) time and 0{n \ogn) 
space, such that for a query point q E P, the largest circle Cq in P, that contains q, can be 
reported in O(logn) time. 
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Proof. The correctness follows from the above discussion. Preprocessing time and space have 
already been established in Lemma [5j We now analyze the query time. The PLiCA query 
requires O(logn) time [16]. If we are in the affirmative case, then finding the node v* at 
the maximum level in T such that q G MEC„* needs another 0(log?T,) time. The QiC query 



for MEC„. can be executed in O(logn) time (Lemma 10). In the negative case, finding the 



appropriate sub- polygon and then performing the QiM query requires O(logn) time. □ 
3.6 Description of the QiC Data Structure 



Recall from Section 3.3 that the QiC data structure preprocesses the medial axis and a 
specified node v such that when queried with a point q G MEC^,, the largest MEC containing q 
can be reported efficiently. 

We use the concept of guiding circles associated with node v E M to find Cq. Let R be an 
array containing the radii of the MECs centered at all the nodes in M, sorted in increasing 
order. 

Definition 1. (Guiding Circles of a node v of M) An MEC C centered somewhere on M is 
called a guiding MEC of the node v of M if (i) its radius is in R, (ii) every MEC on the path 
from V to the center of C in M (both inclusive) is no larger than C, and (Hi) C overlaps 
with MECt,j^ (See Figure^ for an illustration of guiding circles on a single path from v.) We 
denote the set of guiding circles of the node v E M by Sy. 

Computing S„. We can compute S^, by adapting either depth-first search or breadth-first 
search traversal on M starting from v. As we traverse M using (say) depth first search 
starting from v, we keep track of the largest MEC along the path from v to the current 
position in the traversal. When we encounter an MEC C that fits our definition of a guiding 
circle, we include C in along with the id of the mountain in which it is centered. 



Algorithm 4 Preprocessing Phase of QiC 
Require: Polygon P, its medial axis M and a vertex v of M. 
1: Compute the radii of MECs centered at nodes of M and store them in a sorted array R^. 
2: Compute the guiding circles Sy of node v. {We can use an adaptation of either depth 

first search or breadth first search.} 
3: To each C G Sv centered at c, attach the mountain id of the mountain containing c. 
4: For each r G Ry, attach the set SI = {s E Sy\ radius of s is r}. {In Lemmajsj we will see 
that, for any r, \Sy \ is a constant.} 



Before we provide the pseudocode for the query phase, we establish a few lemmas. Recall 
from Lemma [2] that if a guiding circle C contains q, then every guiding circle from MEC^, to 
C will contain q. The proof for the following lemma follows from the definition of guiding 
circles. 

* We say that two circles overlap if they have a common point in their interior. 
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Lemma 7. Let U be the path on M from v to some guiding circle C. 

1. The radii of guiding circles along U are non- decreasing. 

2. Furthermore, ifr&Ris the radius of C and is the radius o/MEC^,, then for every r' G R 
such that Ty < r' < r, there is at least one guiding circle of radius r' in the path from v 
to the center of C . 

Corollary 3. Given a query point q, let pq be the radius in R^, such that 3C G S^^' that 
contains q but Mr > pq, $C G SI that contains q. Then, for every r' G Rt, such that r' < pq, 
3C G that contains q. 



Radii in R are 




X~8jX1S 



" s The arrows heads on 

'Position of v the a;-axis indicate 

centers of guiding circles. 

Fig. 5. For the sake of intuition on tlie construction and usefulness of guiding circles, we show a path 77 laid out on 
the a::-axis. 

Corollary |3] will allow us to perform a binary search for pq which in turn will lead us to the 
largest guiding circles in S„ that contain q. The following lemma ensures that the binary 
search will run in O(logn) time. 

Lemma 8. For any r G Rt,, |SJ^| is bounded by a constant. 

Proof. Consider any r G R. Let be the MECs of radius r in S„. For convenience, let us 
assume that SJ^ does not contain MEC centered at any node of M. Since MECs centered at 
nodes of M have distinct radii, at most one MEC in SJ^ can be centered at a node. 
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(b) 

Fig. 6. (a) Bounding |S|, and (b) Illustration of Lemma|8] 

By the condition (ii) of Definition [T| p?, < r. Also recall that every MEC in Sy must intersect 
MECt, (see Figure [6]^ a)). Therefore, every MEC in must lie entirely within a circle x of radius 
Pv + 2r centered at v. Thus, we need to prove that the number of guiding circles of radius r 
at node v inside x is bounded by a constant. 

Let us consider a point a & P. Let Sl{a) C be the set of MECs in that enclose a. 
Consider any MEC C G S^(a); let c be its center. Let pi and p2 be the two points at which C 
touches the boundary of the polygon P. The chord pip2 must intersect the medial axis (see 
Figure |6|^b)). Note that, the points v and c lie on different sides of piP2- On the contrary, if 

V and c lie in the same side of p'iP2, where p[ and P2 are the points of contact of the said MEC 
and the polygon P, then we can increase the size of the MEC by moving its center c towards 

V along the medial axis (see Figure |6] (b)). Thus, C ^ S^. Thus, we have Zpiap2 > 7r/2. 
Again, the angles subtended by different MECs in SJ^(a) are disjoint. These two facts imply 
that |S^(m)| < 4. In other words, any point inside the circle x can be enclosed by at most 
four different circles of S^. We need to compute \Sl\. Let us consider a function 77(a) defined 
as the number of circles in that overlap at the point a, a G Clearly, rj^a) < 4 for all 
a E X- The total number of circles in can be bounded as follows: 

Total area of circles in ^ f^—zr^. ^)gY vi.^^ — 4vr(r^ -|- 2r)^. 
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Therefore, |s;| < MTe+M! < ^30! = 36. □ 

Answering QiC query Given a query point q and a node f in M such that q G MECt,, we 
compute Cq as follows. Let p G R be the radius of the largest guiding circle in Sy containing q, 
and S(J(g) be all the members of S„ with radius p that contain q. We report Cq by executing 
the steps in Algorithm |5| 



Algorithm 5 Query Phase of QiC 
Require: A query point q lying inside MEC„. 

{We want to find p G R„ such that 3C G that contains q but Vr > p, G that 

contains q.} 

1: Perform a binary search in the array to identify p G Rt,. This also returns the members 
in S(J(g). Note that, each member C G S^(g) is attached with its corresponding mountain- 
id. 

2: For each member in C G S^(g), locate the largest MEC containing q in the mountain 

attached to C by executing the QiM query algorithm. 
3: Report Cq as the largest one among the MECs obtained in Step 2. 



Lemmata P and 10 state the correctness and complexity. 



Lemma 9. At least one of the circles in S(J(g) is centered at some point on the mountain in 
which Cq is centered. 

Proof. Since M'^ is a subtree of M (Lemma |3]), if we explore all the paths in M from node 
V towards its leaves, the center Cq of Cq is reached in one of these paths, say U. Let C be 
the last guiding circle when going from v to Cq. Let the center of C be c'. As a consequence 
of Lemma [7| C G S(^(g). Suppose for the sake of contradiction, C is not centered on the 
same mountain on which Cq is centered. Then, between c' and Cq there is a valley point a on 
the path 77, such that the radius of MEC^ is less than the radius of MECc'. Also, there exists 
another point /3 on the path U between a and Cq such that the radius of MEC^ is equal to the 
radius of MEC„/. Since the radius of MEC^ matches with an element of R, MEC/3 must also be a 
guiding circle. This contradicts our assumption that C is the last guiding circle between v 
and Cq. □ 

Lemma 10. (i) For a node v, can be computed in 0{n\ogn) time and 0{n) space, (ii) 
If the query point q lies in MECi,, then Cq can he computed in 0(log?2) time. 

Proof, (i) First of all, note that |R| G 0{n). The breadth first search in M needs 0{n) time. 
The time for computing the members in S„ is ^^g^ |SJ]| = 0(|R|) (by Lemma|8|, which may 
be 0{n) in the worst case. A sorting of the members in S„ with respect to their radii is 
required; this takes 0(n log n) time. Once sorted, attaching with each r G R will take 
0{n) time. The space requirement can be argued similarly. 
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(ii) If g G MEC^,, the binary search in Sy considers at most 0(log |R|) = 0(log?2) distinct radii. 
For each radius, the number of guiding circles inspected to find whether any one contains 
q is bounded by a constant (see Lemma [s]). Thus p, the largest radius among the guiding 
circles of node v that contains q, can be identified in 0(log?T,) time. 

Let S^(g) denote the set of guiding circles of node v of radius p that contains q. Each of them 
is attached with the corresponding mountain- id. For each member in C G Sf^{q), we invoke 
QiM query to find the largest MEC in the associated mountain M^; this takes 0(log \Mi\) time, 



where may be 0{n) in the worst case (see Subsection 3.2). □ 



Corollary 4. Suppose QiC is restricted to a connected M* C M, i.e., v G M* and M'^ C M* . 
Suppose further that n* edges of P induce the edges in M. Then, the preprocessing time and 
space for QiC are 0(n* log n*) and 0{n*), respectively. The query time will he 0(logn*). 

Proof. We can restrict our R to radii of MECs centered on nodes only in M* . Hence |R| G 0{n*). 
Rest of the proof follows from the previous discussion. □ 



Proof (of Theorem^. The proof follows from Lemma [o] and Lemma 10 



□ 



4 QMEC problem for Point Set 

The input consists of a set of points P = {pi,P2, ■ ■ ■ , Pn} in I^^- The objective is to preprocess 
P such that given any arbitrary query point g G M^, the largest circle Cq that does not contain 
any point of P but contains q, can be reported quickly. Observe that, if q does not lie in the 
interior of the convex hull of P, then we can easily report a circle of infinite radius passing 
through q, that does not overlap with P. So, in the rest of this section, we shall consider the 
case where q lies in the interior of the convex hull of P. 

Consider the Voronoi diagram of P. Observe that the MEC centered at any Voronoi vertex 
touches at least three points of P. To simplify our presentation, we assume that MECs centered 
at Voronoi vertices are of distinct sizes. In the course of our algorithm, we treat the Voronoi 
diagram of P, as a plane graph G. To keep G within a finite region, we insert artificial 
vertices, one for each unbounded edge in the Voronoi diagram of P, so that G is the plane 
graph of the Voronoi diagram of P with each unbounded edge clipped at its corresponding 
artificial vertex. In placing the artificial vertices, we ensure that (i) every MEC centered at an 
artificial vertex must be larger than all the MECs centered at Voronoi vertices, and (ii) the 
MECs centered at artificial vertices do not overlap pairwise within the convex hull of P. They 
may overlap outside the convex hull of P. The second condition ensures that there exists no 
query point q, in the convex hull of P, which can be enclosed by more than one MEC centered 
at artificial vertices. From now onwards, we will use the term vertices of G to collectively 
refer to Voronoi vertices and artificial vertices. We will use both the geometric and graph 
theoretic properties of G. In particular, to achieve the subquadratic preprocessing time, we 
use the classical planar separator theorem [22] . The intuition is as follows. 
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Consider the following naive approach to solving QMEC on points. Suppose we store the MECs 
of vertices in G in a PLiCA data structure. Suppose, furthermore, that we preprocess each 
vertex for QiC adapted for points set. We note that here also the QiC data structure of 



a vertex can be implemented using guiding circles (cf. Section 4.1 for details). Given an 
arbitrary query point q in the convex hull of P, we know that q lies in at least one of the 
MECs centered on vertice^ so we can locate one such MEC, say MEC^, for some vertex v. We 
can execute the QiC query for the point set to identify Cq. This, unfortunately, will require 
0(n^ log n) time for preprocessing because each QiC preprocessing requires O(nlogn) time. 
Instead, to achieve sub-quadratic bounds on the preprocessing time and space for the QMEC 
problem we employ a divide-and-conquer approach by recursively splitting the vertices of G 
using the planar separator theorem stated below. 

Theorem 4. A planar graph G on n vertices can, in 0{n) time, he partitioned into 
disjoint vertex sets A, B, and W such that (i) \W\ G 0{y/n), (ii) \A\, \B\ < 2n/3, and (Hi) 
there is no edge in G that joins a vertex in A to a vertex in B. 



Algorithm 6 Preprocessing Phase of QMEC for set P of points in 

Require: This is a recursive algorithm. In the first call, the input graph Gin is G. Subse- 
quently, the input graph Gin is a subgraph of G. 
Ensure: The first call on the entire plane graph G will return a pointer r to the root of the 
separator decomposition tree T. 
1: if Gin is empty then 
2: return NULL. 
3: end if 

4: Create a node v of the separator decomposition tree T. 

5: Compute the planar separator vertices Wi^ of Gin (cf. Theorem |4]). Denote the two 
separated subgraphs as A and B. 

6: Compute a PLiCA data structure <P on the MECs centered on vertices in Win_- 

7: Compute a PLiCA data structure on the MECs centered on vertices in A. 

8: Compute a QiC data structure corresponding to node v. {This QiC data structure is built 
on the two promises that (i) at least one MEC centered on the planar separator vertices 
Win will enclose the query point g, and (ii) Cq is centered on an edge of the plane graph 

Gin.} 

9: Attach ^, and the QiC data structures to v. 
10: f.LEFxCHlLD ^ (Call Algorithm [6] on A). 
11: ?;.RightChild ^ (Call Algorithmje] on B). 
12: return pointer to v. 



The MECs on vertices can be thought of as circumcircles of triangles in the Delaunay triangulation of P and therefore 
the union of these MECs covers the entire convex hull region. 
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Algorithm 7 Query Phase of QMEC for set P of points in 

Require: A query point q inside the convex hull of P and pointer r to root of the separator 

decomposition tree T. 
Ensure: The larges MEC Cq that contains q is returned. 
1: PTR <— r. 

2: while PTR is not null do 

3: Let t G T be the node pointed by PTR. 

4: Let Gt be the subgraph of G associated with t. 

5: Let Wt be the separator vertices of G. 

6: Let At and Bt be the two disconnected subgraphs obtained when vertices in Wt are 
removed from Gt- 

7: Let $t, 6*4, and QiCj be the data structures attached to t. 

8: if 3w e Wt such that w encloses q (we check this using <Pt) then 

9: Gq -ir- (circlc rctumcd by querying QiC^ with query point q). 
10: return Gq. 
11: end if 

12: if 3w' in the data structure associated with v such that w' encloses q then 

13: PTR = i;.LEFTCHILD. 

14: else 

15: PTR = W.RIGHTCHILD. 

16: end if 
17: end while 

18: {The execution will not reach this point.} 



We construct a separator decomposition tree T as follows (cf. Algorithm |6] for a detailed 
pseudocode). The root r of T represents the plane graph G. We attach two PLiCA data 
structures <P and O at r. In we store MECs centered on the 0{^/n) planar separator 
vertices (denoted by W). We also build the QiC data structure for the node r. The details of 
QiC in the current context of a set of points (rather than polygon) is described in Subsection 



4.1[ For now, however, we state the QiC problem in the current context where P is a set 
of points. The node r has 0{^/n) MECs corresponding to the 0{\/n) separator vertices and 
therefore, the QiC attached to r comes with two promises. The first promise is that the query 
point q will be enclosed by at least one of the separator MECs. (Note that this first promise 
is an adaptation from the context where P is a simple polygon. In that context, because the 
medial axis of P was a tree, the separator was a single vertex.) Our second promise is that 
Gq is centered on some edge of the plane graph G attached to r. In the query phase of QiC, 
given a query point g, we are to return Gq. 

The removal of the vertices in W from G will induce two disjoint subgraphs A and B. Without 
loss of generality, we pick A and build containing the MECs centered at the vertices of A. 
The root r has two children LeftChild and RightChild in T. LeftChild is associated 
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with the subgraph A while RightChild is associated with the subgraph B. The two children 
of V are then processed recursively. 

In the query phase (cf. Algorithm [?]) , we are given a query point q. We find the highest node 
t in T such that (at least) one of the MECs C stored in the associated ^ encloses q. We find t 
by a traversal from the root node of T. Let v be the center of C. The point w is a separator 
vertex in the graph Gt associated with t. Recall that each separator vertex has a QiC data 
structure (restricted to Gt) associated with it. We prove subsequently that when we query 
the QiC data structure attached to v with the query point q, we will indeed obtain the largest 
MEC Gq that encloses q. 

We now turn our attention to analyzing Algorithm |6] and Algorithm [7} We begin with some 
important lemmata. 

Lemma 11. Consider any cycle H in the Voronoi diagram of P. Let Gh be an MEC centered 
at some point on H . Then, there exists another MEC G'jj centered at some other point on H 
that does not properly overlap with Gh- 

Proof. Clearly, any cycle in the Voronoi diagram of P must contain at least one point from P 
inside it. Let p G P be such a point that lies inside the cycle H (see Figure [?]). Let Gh be any 
MEC centered at some point on H; let ch be the center of Gh- Consider the line connecting 
ch and p. It intersects H at another point c'^. It is easy to see that the MEC C^, centered at 
d^, will not properly overlap Gh as, otherwise, p will he inside both Gh and G'^. 




Fig. 7. Illustration for Lemma [TT] 
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Lemma 12. (Unique Path Lemma.) If C and C are two distinct but overlapping MECs with 
centers at c and d , respectively, then there is a unique path n{c, c') from c to d in the Voronoi 
diagram of P such that every MEC centered on that path encloses C (1 C . 



Proof. The structure of the proof is as follows. We provide a procedure that constructs a 
path n{c, d) from c to d along the Voronoi edges, and ensure that every MEC centered on that 
path encloses CdC. As a consequence of Lemma 11 , the path does not form an intermediate 
cycle and terminates at d. Finally, we again use Lemma 11 to show that no path 77', other 
than n{c,d), exists between c and d such that every MEC centered on U' contains C fl C. 
Throughout this proof, we closely follow Figure |8] in order to keep the arguments intuitive. 
To keep arguments simple, we assume that c and d are Voronoi vertices. The arguments hold 
even when c and d are not Voronoi vertices. 

Let a be the number of points in P that C touches. These a points partition C into a 
arcs. The degree of the corresponding Voronoi vertex c (center of C) is also a because each 
adjacent pair of points of P on the boundary of C will induce a Voronoi edge incident on c 
and vice verse. These Voronoi edges and their corresponding arcs are denoted by Cq and Sq, 
for 1 < i < a. 

Consider the other MEC C C and centered at a vertex d) that overlaps with C. C 
intersects C at two points ti and ^2- Since C is empty, both ti and t2 must lie on one of the 
a arcs of C . Let us name this arc by s;^. Consider the edge = (c, C2) that corresponds to 
the arc Sq. The other end of e^, i.e., the vertex C2, is called the next step from c toward d 



and denote it as ns(c, c'). Consider the pseudocode in Procedure |8] that generates the path 
denoted by i7(c, d): 



Algorithm 8 i7(c, d) Computation 

1: n{c, d) ^ (c) 

2: next ^ c 

3: repeat 

4: next ^ ns(next, d) 
5: Append next to i7(c, c') 

6: until next equals d. {Note that this is the only terminating condition.} 



We now show that (i) i7(c, d) is our desired path, and (ii) there exists no other path satisfying 
the unique path lemma. 

Proof of correctness: Algorithmjsjconstructs a path i7(c, d) = (ci = c, C2, . . . , Q, Q+i, . . . ,d), 
where each Cj is a vertex in the Voronoi diagram of P. Let C2 denote the MEC centered at C2. 
If C2 = C, then the procedure terminates and, as required, every MEC centered on the edge 
(c, C2) encloses C n C2 = C n C". 

Therefore, we consider the case where C2 7^ C". We need to prove C fl C" C C fl C2. 
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Let pi,P2 G P be the points at which C and C2 intersect; pi,P2 are the end points of the 
arc s;^ that defines the move of the next step toward c' (in Figure [sj j is 2). By definition, 
ti and ^2 he on the arc s;^. Notice that C H C" (shaded region in Figure [s]) is shaped hke a 
rugby ball with ti and ^2 at its end-points. One side of C H C (called the initial side) is in 
C and the other side (called the final side) is in C. Clearly, ti and ^2 are inside (or on the 
boundary of) every MEC centered on the edge e^. Otherwise, as we go from C to C2, a circle 
would be present that must touch the final side of CCiC, but that would mean that we have 
either 

— reached C, which contradicts our assumption that C2 ^ C, 

— or found a MEC that contains C, which contradicts the fact that C is itself a MEC. 

We now make two observations: (i) C touches the initial side of C fl C, but (ii) no other MEC 
centered on e-jj (C2 in particular) touches the final side of C H C. 

Observation (i) is obvious. We prove Observation (ii) by contradiction. Let C* be an MEC 
centered on e-^ that touches the final side of C H C at, say, some point t*. It is easy to see 
that C* will contain C because C* touches C at the point t* and also C* contains ti and 
t2, which are also on the boundary of C . Thus we have a contradiction that C is an MEC. 
Thus, it is clear that C n C" C C n C2. 

Consider two adjacent vertices c, and Cj+i along i7(c, c') with MECs Cj and Cj+i centered on 
them, respectively. The above argument can be easily extended to give us the following: 

anc c a+i n c. 



Therefore, we can conclude that every MEC along i7(c, c') encloses CdC. Lemma 11 suggests 
that n{c,c') does not form a cycle. The only stopping condition is when we actually reach 
c'. So n{c, c') terminates at c' in at most 0{n) steps. 

Proof of uniqueness: To complete the proof of this lemma, we must show that i7(c, c') is 
the only required path. For the sake of contradiction, assume that there is another path U' 
such that every MEC centered on U' contains C (iC Then, there are two distinct paths from 
c to c' such that every MEC centered on both the paths contain C fl C. Clearly, there must 



be a cycle when the two paths are combined. From Lemma 11, we know that there are pairs 
of MECs in the cycle that do not overlap on each other. This is a contradiction. Thus i7(c, c') 
is the only required path. □ 

Recall that, given a query point q, we locate Cg by traversing the tree T from its root node 
r. At each node t on the search path, we search in the $t data structure to check whether 
q lies in an MECs corresponding to a separator vertex of node t. If there exists an MEC„ G $t 
containing q, then we perform QiC query in t to identify Cg. Otherwise, we search q in 0t, 
associated with the partition Af. Now, if there exists an MEC„ G Ot containing q, we proceed 
towards the left child of t, otherwise we proceed towards the right child. 

Lemma 13. The search with q must stop at a node t of T, and outputs a vertex v in the 
plane graph Gt associated with t such that (i) q & MEC„ and (ii) Cq is centered on some edge 
ofGt. 
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Proof. Because MECs on Voronoi vertices are circumcircles of triangles in the Delaunay trian- 
gulation of P, the union of these MECs covers the entire convex hull region. Since every MEC 
centered on a Voronoi vertex of G is a separator for some node in T, the proof of (i) follows. 

Suppose some MEC C centered somewhere outside Gt encloses q. From the planar decompo- 
sition of G down to Gt, it is clear that the (collective) neighborhood r{Gt) of Gt consists of 
vertices that appear in the separator vertices associated with some ancestor of t. Therefore, 
by the Unique Path Lemma, there must exist a vertex v* G r{Gt) such that the MEC centered 
on V* encloses q. Since v* is associated with some ancestor t* of t, the search path in T must 
have stopped at t* instead of coming all the way to t, which establishes a contradiction. □ 



Now that Lemma 13 is established, we can easily see that if the search path in T for a query 



point q stops at node t G T, then the two promises required for QiC^ data structure are 



fulfilled. Therefore, assuming QiC is correctly designed (established in Section 4.1), we get 
the following lemma. 

Lemma 14. Algorithm preprocesses a set P of points such that given a query point q, 
Algorithm^ can he used to correctly finds Cg. 

In the next subsection we describe both the preprocessing and query phases of QiC before 



proving time and space bounds in subsection 4.2 



4.1 QiC data structure for points set 

The QiC data structure for the points set case (attached to nodes in T) largely mimics the 
simple polygon case. Note that each t G T has a QiC data structure attached to it. We 
reiterate that, while the query point q is promised to lie in a particular MEC in the polygon 
case, in the points set case, q is promised to lie in at least one of the MECs centered on a node 
in Wt- The preprocessing and query algorithms are given as self explanatory pseudocode in 



AlgorithmPand Algorithm 10, respectively. The time and space bounds of the preprocessing 



phase follow from the following lemma. 

Lemma 15. For any node t E T, any vertex v in Gt and any r E R, we define = {C \ 
C E A and radius of C is r}. We claim that \Sl\ is bounded by a constant. 

Proof. The key ideas required to prove this lemma have already been discussed in the context 
of Lemma |8| Therefore, we limit ourselves to making a few important observations that 
establish a correspondence between the current context (where P is a set of points) to the 
context of Lemma [s] (where P is a simple polygon). 

Firstly, observe that all circles in must lie in a circle x of radius + 2r centered at v] 
here is the radius of MEC^ and r is the radius of circles in (see Figure [6]). 

To make the second observation, consider a circle C E centered at a point c strictly in 
the interior of an edge e = (fi,f2). Without loss of generality, assume MEC^,^ is smaller than 
MECt,2. Therefore, C will be no smaller than MEC„j and no larger than MEC^j. Our second 
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Algorithm 9 Preprocessing for QiC attached to node t G T. 

Require: A node t G T and the plane graph Gf attached to t along with separator vertices 
Wt and parts At and Bf. 
1: Compute R = {r I 3 an MEC of radius r centered on a vertex of Gt}. 
2: for all ^; G Wt do 

3: Compute S^, = {C I C is an MEC that fulfills the following}: 

1. C overlaps with MEC„, 

2. Radius of C is in R, and 

3. Every MEC in the unique path from MEC^ to C has radius no more than that of 
C. 

4: Sort 

5: Compute BFS tree A^, rooted at v with centers of MECs in S^, as nodes. 

6: For each node u in A^,, there is at most one edge r] incident on v such that the MECs 

centered on rj strictly grow in size (starting from i/). Mark rj red. 
7: end for 

8: The QiC data structure associated with t consists of and Ay for all v G Wt- 



Algorithm 10 Query phase for QiC attached to node t G T. 
Require: QiC data structure and a query point q that meets the two promises. 
Ensure: The largest MEC containing q is computed and returned. 
1: Using the PLiCA data structure ^ attached to t, we find a,v eWt such that MEC^ encloses 
Q- 

2: Perform a binary search on to find the radius r^ax of the largest MEC in S„ that encloses 
q. Also compute 

C = {C I (C G S^) A (radius of C = rmax) A (C encloses q)}. 

3: Cmax C, where C G C is chosen arbitrarily. 
4: for all C G C do 

5: Let c be the center of C. Recall that c is a node in 

6: Let C* be the largest MEC centered on the red edge in Ay incident on c. If no red edge 

is incident on c, assign C* <^ C. 
7: if C* is larger than Cmax then 

9: end if 
10: end for 

11: return Cmax- 



observation is that the MECs centered on e are growing in size in the vicinity of c as we move 
in the direction from Vi to V2- 
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To make the third observation, note first that the unique path (as defined in the Unique 
Path Lemma) from v to c passes through vi and not through V2- Let pi and p2 be the two 
points in P that touch C. The third observation is that the chord pip2 intersects the unique 
path (as defined in the Unique Path Lemma) between v and c (see Figure [6] for a similar 
situation in the context where P is a simple polygon). 

The rest of the proof follows from the proof of Lemma |8] □ 

Lemma 16. Algorithms\^ and [7^ correctly implement the QiC data structure. The prepro- 
cessing time and space of the QiC attached to node t in T is bounded by O(n^^^lognt) and 
0{n^J'^), respectively, where rit is the number of vertices in the plane graph Gt attached to t. 
Queries can be answered correctly in O(lognt) time. 

Proof. Let Gt be the subgraph attached to a node t G T, and Wt be the separator vertices of 
Gt- Recall that \Wt\ = 0{^^/\r^~\), where nj = \Gt\- In the QiC data structure for the node t, 
\Sy \ can be 0{nt) for each node in f G Wt, and it can be computed in 0(nj logn^) time. Thus, 
the time required to create the QiC data structure for all the nodes in Wt is 0{n^^'^ log Ut). 
The space requirement is 0{n^^'^). 

The correctness argument is similar to the polygon case. In the polygon case, a single path 
between any two points on the medial axis followed immediately from the fact that the 
medial axis was a tree. In the current context, the Unique Path Lemma provides a similar 
unique path. □ 



4.2 Complexity 

Lemma[T4ljustifies that our proposed algorithm correctly computes the largest MEC containing 
the query point q among the points in P. The following lemma establishes the complexity. 

Lemma 17. The preprocessing time and space complexities of the QMEC problem are 0{n^ log^ n) 
and 0(n2 logn), respectively. The query can be answered in O(log^n) time. 

Proof. The preprocessing consists of the following steps: 

• Constructing the tree T. At each node t of T, we need to compute the separator vertices 
among the set of vertices of the Voronoi subgraph Gt corresponding to the node t. The 
time complexity for this computation is 0(|Vt|), where \ Vt\ denotes the number of vertices 
in Gt. Since the total number of vertices at each level of T is 0{n), the total time spent 
for computing the separator vertices at all nodes in each level of T is 0{n). Since the 
height of T is at most O(logn), the total time for constructing it is O(nlogn). 

• For each node t in T, we need to construct the QiC data structure. The subgraphs associ- 
ated with each node in any particular level of T are disjoint. Therefore, as a consequence 



of Lemma 16, the preprocessing time and space required for QiC data structures associ- 
ated with nodes in any particular level is 0(n^/^logn) and 0{n^^'^), respectively. Since T 
can have at most O(logn) levels, the preprocessing time and space complexities follow . 
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• For each node t, we also will need to compute two PLiCA data structures <P and O, but 
the time and space complexities of the QiC data structures dominate the complexities of 
computing the PLiCA data structures. 

While querying with a point q, searching in the PLiCA data structures ^ for each node in 
the search path of T will take O(logn) time. We have to traverse a path of length at most 
O(logn) to get to a node t in T such that there is a vertex v of Gt such that q e MEC„. Thus, 
traversing T needs O(log^n) time. Finally searching in Sy and to get Cg needs another 
O(logn) time (see Theorem [3]). □ 

4.3 Improving the query time 

We now show that a minor tailoring of the data structure reduces the query time to 
O(lognloglogn), while maintaining the same preprocessing time and space. 

Data structure. After computing the planar separator tree T, each MEC C centered on a 
vertex in G is attached with 

— an id, which is the level of T in which C belongs as a separator MEC, and 

— a pointer to the node t in T such that C belongs to the separator vertices of Gt- 

Next, we create an array F of O(logn) data structures as follows. Each Fi is a PLiCA data 
structure constructed with the set of MECs with id ranging from 1 to i, i.e., root to the level 
i. 

Query. While querying with a point q, we conduct a binary search on the array F of data 
structures to find F^ such that there is an MEC C* in F^ that contains q, but no MEC in Fi_i 
that contains q. Let t* the node in T where the center of C* is a separator vertex. We now 
perform a QiC query on t* and report the result of that QiC query as the required MEC Cg. 

Theorem 5. The improvement described in this subsection is correct and its preprocess- 
ing time and space complexities for the QMEC problem are 0{n^^^\og^ n) and 0{n^^'^\ogn) 
respectively. Each query can be answered in O (log n log log n) time. 

Proof. The correctness follows from the fact that there is no MEC with id smaller than C* 
that contains q, but C* in fact contains q. Therefore, a QiC on C* indeed gives us the required 

MEC Cg. 

Each Fi requires 0{nlogn) time and 0{n) space. Therefore, to construct F, we require 
0(n log n) time and 0(n log n) space, which are subsumed in the bounds established in 
Lemma [17] to construct T. 

In the query phase, each PLiCA query on any element of F requires 0(log?7.) time and 
the binary search over all elements of F requires O(loglogn) such PLiCA queries, thereby 
requiring O(logrzloglogn) time overall. The QiC query requires an additional O(logn) time, 
which is subsumed. 
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4.4 Achieving O(logn) Query Time 

Here, we shall use Frederickson's r-partitioning of planar graphs, stated below, to improve 
the query time complexity to O(logn). Furthermore, this algorithm is simpler in that it does 
not require us to construct a divide and conquer tree. 

Lemma 18. [T4^ Given a planar graph G with n vertices with a planar embedding and a 
parameter r (1 < r < n), 

(a) G can he partitioned into 0{-) parts with at most 0{r) vertices in each part, and a total 
of 0{^) boundary vertices over all the partitions. 

(b) This partitioning can be computed in 0[nlogn) time. 

We compute the r-partitioning of the graph G with r set to n^/^. Now, we construct two 
data structures, T and \P, as stated below. 

T: We construct a PLiCA data structure and a QiC data structure over all MECs centered on 

boundary vertices in the r-partitioning. 
^: It consists of a PLiCA data structures with the set of MECs that correspond to the internal 

vertices of all the partitions. Furthermore, for each partition j, we construct a QiC data 

structure limited to partition j on the MECs centered on the internal vertices of partition 

J- 

For a given query point q, we first search in T to check whether there exists a boundary MEC 
that contains q. Here two cases may arise: 

Case 1: If an MEC C G T encloses q, then we search in the QiC data structure attached with 

T to identify the largest MEC containing q. 
Case 2: Otherwise, we search in ^ to identify an MEC G' that contains q. We also find the 

partition j on which G' is centered. We then search in the QiC data structure attached 

with j' to identify the largest MEC containing q. 

Lemma 19. The above algorithm correctly identified the largest MEC containing q. The pre- 
processing time, space and query time complexities of this algorithm are 0(n^^^ logn), 0{n^^^) 
and 0(\ogn) , respectively. 

Proof. The correctness of the algorithms follows from the following argument. During the 
query, if Case (i) arises the algorithm produces the correct result since the QiC data structure 
attached to T is built on the MECs corresponding to all the vertices in G. If Case (ii) arises, 
and q lies in an MEC G' of the j-th partition, then it implies that q lies in some MEC in the 
proper interior of the j-th partition. Thus, the largest MEC containing q is surely an MEC G* 
of the j-th partition. Thus the QiC of G' constructed with the MECs in the j-th partition only 
is sufficient to obtain Gg. 

Now, we justify the complexity results of the algorithm. The total size of the QiC data struc- 

2 2 

tures in T is O(^), and these are constructed in O(^logn) time. The total size of the QiC 



27 



data structures for all the MECs in the j-th partition of \I/ is O(r^), and these are constructed 
in 0(r^ log r) time. Since, we have at most O(^) partitions, the total space and time required 
to construct ^ is 0{nr) and (nrlogr), respectively. Thus, the total preprocessing space and 
time complexities are 0{^ + nr) and 0{{^+nr) logn), respectively. Choosing r = 0{'n?/^)^ 
the preprocessing time and space complexity results follow. 

The query time complexity follows from the fact that the search in the PLiCA of both T 
and ^ take 0(log?T,) time, and the search in the QiC data structure of exactly one MEC needs 
another O(logn) time in the worst case. □ 
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Fig. 8. Illustration of i7(c, c') in unique path lemma 
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